Spin Density Waves and Domain Wall Interactions in Nanowires 
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We investigated how the dynamics of a domain wall are affected by the presence of spin density 
waves in a ferromagnetic wire. Domain walls and other scattering centres can cause coherent spin 
density waves to propagate through a wire when a current is applied. In some cases the spin torque 
due to these scattered electrons can be enhanced such that it is on a par with the exchange and 
anisotropy energies controlling the shape and dynamics of the domain wall. In such a case we find 
that the spin density waves enhance the current induced domain wall motion, allowing for domain 
wall motion with smaller current pulses. Here we consider a system involving two domain walls and 
focus on how the motion of the second domain wall is modified by the spin density waves caused by 
the presence of the first domain wall. 

PACS numbers: 75.60.Ch, 75.60.Jk, 75.75.+a, 75.30.Fv 



I. INTRODUCTION 

The dynamics of a domain wall (DW) in a ferromag- 
netic wire has received much interest due to both its fun- 
damental physical importance, but also in light of pos- 
sible future applications^— Coupling between the carri- 
ers of spin polarized currents and the magnetic moments 
forming DW alters not only the transport properties^— 
of the wire but also the magnetization itself. The cur- 
rent traveling through the wire couples to the domain 
wall causing it to move through the wire. Thus, differ- 
ently orientated collinear regions of the wire may serve 
as discrete bits that can be steered by the current. In 
order to increase the efficiency of these devices^ it is 
necessary to increase the density of the DWs. Therefore 
even a weak interaction between DWs can become im- 
portant. This coupling can also distort the DW during 
its propagation through the wire. Some examples from 
the now voluminous literature can be found herei 4 ' 9 ' 10 

So far most work has concentrated on how the spin 
polarized currents will affect the domain walls. How- 
ever, there are also feedback processes which can lead to 
interactions between the DWs and spin density fluctu- 
ations. The relevant effects are enhanced considerably 
when considering low dimensional structures*^ An effec- 
tive RKKY interaction between two domain walls in a 
wire can be mediated by the current electrons J2r— The 
energy profile of this interaction then favours particu- 
lar alignments and positions of the domain walls. Fur- 
thermore any change in the spin density caused by the 
presence of one domain wall can have an effect on a wall 
further down the wire, defined in the direction of electron 
flow. 

For applications it is most interesting to consider a 
high density of DWs in a wire. In this case it is crucial 



to inspect how the DWs affect each other when a spin 
polarized current is sent down the wire. In a previous 
article^ we looked at how the current mediated RKKY- 
like interaction changes the magnetization dynamics with 
a focus on the motion of relatively sharp walls. These 
walls were treated as localized moments. Here we extend 
this idea to consider how a long (adiabatic) domain wall 
is distorted in the presence of a spin density wave caused 
by scattering from a previous domain wall. The spin 
density wave acts as an effective applied magnetic field. 
However, because the spin torque of the current electrons 
is no longer uniform, it is not a homogeneous field. This 
will tend to distort the shape of the domain walls. 

Furthermore this non-uniform spin density allows the 
domain wall to be set into motion using a smaller current 
pulse than would otherwise be the case. This is closely 
related to its effect as an effective magnetic field. This of- 
fers an alternative possible solution-^ to the problems of 
overheating associated with sending too large a current 
through the nanowires. Though we note that in order to 
obtain a spin density of the appropriate order of magni- 
tude one must consider very narrow wires, of a few atoms 
across. 

It is now possible to fabricate ferromagnetic wires com- 
posed of single atomic sites on a latticeJ£~— Surpris- 
ingly these wires show both ferromagnetic order and con- 
tain regions of non-collinearity^ in other words: domain 
walls. In the low dimensional limit of these metallic wires 
the effect of the spin density corrections on the domain 
wall may become on the order of that of the exchange 
and magnetic anisotropy. This can cause severe distor- 
tion of the domain wall profile and alter the way it travels 
through the wire. We note that the spin density correc- 
tions must not of necessity originate from a DW, any 
spin coherent scattering centre will lead to the behaviour 
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described in this article, for example a magnetic impu- 
rity or region of non-collinearity in an otherwise single 
domain ferromagnct. Although here we concentrate on 
the case of a sharp domain wall as the scatterer, the 
generalization to other scenarios is straightforward. The 
spin density for a specific material can be varied by two 
parameters, the overall and relative magnitudes of the 
two spin channels. By varying these parameters we can 
consider a variety of scattcrcrs. 

Our set up thus consists of a ferromagnetic wire with 
two domain walls present. A current is then applied to 
the wire. We focus on how the motion of the second 
domain wall (defined with respect to the current direc- 
tion) is affected by the spin density waves caused by the 
presence of the first domain wall in the wire. 



II. THEORETICAL FORMULATION 

We start with a zero temperature general Hamilto- 
nian for a nanowire describing non-interacting conduc- 
tion electrons of spin a, aj,(r), coupled with a strength 
J to some non- uniform bulk magnetization, M (r) i 21 i 22 



where 



H' = J dr^4(r)[^a/3 - Ja c 



p.M{r)}ap{v). (1) 



£ is the kinetic energy operator. The inhomogeneity 
can be partially dealt with by a local vector gauge 
transformation!^— After this transformation we will 
have a uniform Zeeman splitting term and a spin- 
dependent spatially varying potential, U a p{r), which de- 
scribes the scattering from any non-collinear configu- 
rations of the magnetization. This transformation de- 
scribes a local rotation in spin-space to align the magne- 
tization direction throughout the wire. It is possible for 
any form of M(r) which has a constant magnitude. The 
gauge transformation is 
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defined such that (r) er ■ n(r) T(r) = a z , where n is the 
unit vector along Our new Hamiltonian is then 

H = J dr^fit (r)[&,0 - JMoZp + U a0 (r)]a p (r)(3) 

with the scattering potential given by 

U(r) = -^-[2A(r).d r + d r .A + A 2 (r)], (4) 

and A(r) = T^(r) V r T(r) is a gauge potential. 

For the case of two domain walls in a nanowire with 
the first DW located at z = and the second DW at 
z = zq, and of widths L' and L respectively, we have 



^cos[8(z)} sin[(/?(z)] N 
M{z) = M | sm[6(z)} sin[v?(z)] 
cos[<^(z)] 



(5) 



(p(z) = 7r — cos 1 [ tanhfz/Z/]] 
f 7r — cos -1 [tanh[(z — z a )/L]\ . 



(6) 
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The angle 9(z) is used to give the walls a different ori- 
entation. Around the first wall wc set it arbitrarily to 
and around the second to an angle 8 , which therefore 
defines the relative orientation of the two walls. Cross 
terms between the walls are very small and can be ne- 
glected provided Zq L, L' , allowing us to approximate 
the scattering potential, 



V(z) »Ui(z)+U 2 (4 
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as independent contributions from each DW. To see some 
interaction effect between the DWs we must nonetheless 
of course have the distance between the DWs as less than 
the spin coherence lengthscale in the system. Now 
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If L > Xp then wc can also neglect the y 2 ' ( z ) an d [v 3 ^ 2 )] 2 
terms which arc of order of (Xf/L) 2 . 

The correction to the spin-density of a single scattered 
wave of spin 6 is 



AS S = 6^9^. 



(10) 



8tp e s is the scattered wavefunction, from an incoming 
particle of spin 8, at an energy e. It can be calculated 
within the Born approximation provided that V > Xf ^ 
This tells us the nature of the spin waves in our system. 
These spin density corrections result in an inter-DW spin 
torque and RKKY-like interaction, details of which can 
be found in Sedlmayr et alJ^~— . The result for the spin 
density around the second domain wall at an arbitrar- 
ily chosen configuration, 9q = 7r/4, is shown in figure 
H Plotted for A F = 0.367nm, L' = Xp, L = 10X F , 
and JM = 2.24eV. The adiabatic change in spin density 
across the width of the second DW is clearly visible, and 
is modified by the presence of the spin density waves. 
These oscillations in the spin density will tend to inho- 
mogeneously distort the domain wall profile away from 
its equilibrium profile. 



III. MOTION OF THE SECOND DOMAIN 
WALL 

We wish to find the effect of the scattered electrons 
on the second domain wall. We take the spin density 



3 



5- 



c/T 



— S 

__ s 




-2.5 



2.5 



(z-z Q )/L 



FIG. 1. (Color online) The spin density of the carrier electrons 
around the second domain wall, centred at Zo and measured 
in carrier spin per unit length. One can clearly see both the 
spin density waves and the overall profile of the second domain 
wall. The solid (black) curve is for S x (z), the dashed (red) is 
for S y (z), and the dotted (blue) is for S z (z). 



of the scattered electrons and calculate its effect on the 
second DW using the now standard Landau-Lifschitz- 
Gilbcrt (LLG) equation approach.! 21 ' 27 ' 28 In our case this 
gives us 



d t M{z,t) = - 



1 



I + a 2 



jM{z,t) x n 



+ J L M (z,t) x (M(z,t)xH) 



M 2 



M(z,t) x (M(z,t) x (J.V)MOM)) 



bj ^ M a) M(z,t)x(j.V)M(z,t) 



(11) 



7 is the gyromagnetic ratio: 7 = ^j^-, where g is the 
Lande factor, and a is a phenomenological constant char- 
acterizing the magnetic damping of the system. The con- 
stants for the non-adiabatic terms are: j is the current 
density; £ = T ex /r s f is the ratio between the exchange 
and spin-flip time scales; and bj = ^nTj^ffpJ' w ith P the 
polarization. The effective field is 



H(z) = JS(z) + a c 



d 2 M{z) 
dz 2 



+ KM x (z)x. (12) 



a cx is the exchange coupling. The coefficient K char- 
acterizes the magnetic anisotropy in the wire. We have 
chosen here x as the anisotropy axis. We us o 21 ' 22 a cx = 
J /Ma and L2 = it \f K / a eyi . (a is the lattice spacing.) 
Note however that the domain wall width can be changed 
by geometric effects as well as by changing the anisotropy 
or the exchange energy. The effective field H, which the 
DW experiences is composed of the spin density fluctua- 
tions of the carriers, the exchange interaction of the DW 
itself, and the magnetic anisotropy of the system. The 
exchange interaction tends to oppose distortions on the 



wire and will impose an upper limit on how much the 
DW can be distorted by a given spin density. 

We calculate the dynamics resulting from equation 
([TT]) using the boundary conditions M(0, z) = M n (z) and 
M(t, zq ± zi,) = Mq(z ± Zb), where Mq(z) is the magne- 
tization of equation ([5]). z^ must be taken to be large 
enough to insert no artifacts in the results. Note that 
from equation (fT2"j) it is clear that for a sizeable effect 



we require aL > \S\, where S is the spin density. This 
can be reached by considering the experimentally real- 
izable limit of a chain of atoms with a cross section of 
approximately a cs — \ 2 F . The lattice spacing is taken as 
a = 1 nm and we use here the Fermi wavelength of iron: 
\p = 0.367 nm. Also we take M = 1.72 x 10 6 Am" 1 , 
JM = 2.24 eV,je = -2.33 xlO 8 Am ~ 2 , and £ = O.Oll^ 8 . 
The second domain wall has a width of L = 10 Xp, the 
first of the order of the Fermi wavelength, V = Xf- The 
DW we consider the motion of is therefore relatively nar- 
row. This gives us the advantage of making our results 
clearer. However it could also limit the applicability of 
the model. We note though that a DW would be nar- 
rower in a wire with a small cross section than in the 
bulk or thin film case, meaning that upon patterning of 
the wire the DW width can be decreased. Nonetheless 
the results here would hold for DW widths of the order of 
lOnm. The case for DWs of larger widths is discussed be- 
low. Finally: a = 0.01, and the exchange and anisotropy 
are as defined above (x is therefore a hard magnetic axis). 
The second DWs configuration is begun at an arbitrary 
angle, #0 = 7r /4- Altering the initial angle of the second 
DW does not significantly affect the current induced mo- 
tion of the second domain wall. The speed of the DWs 
motion is the same irrespective of the value of 9q. Qual- 
itatively the same motion is seen except for some small 
distortions in the DW profile which depend on 0q, 

Firstly let us consider the effects of the spin density 
wave without considering the current induced motion in 
equation (jlip , plotted in figures [2] and [3] In the right 
hand side (RHS) of figures [5] and [3] over a short period of 
time the wall is distorted into a low energy position, as 
in the case without the spin density corrections (the left 
hand side (LHS) of figures [2] and [3]) , but as the spin den- 
sity fluctuates on a scale much shorter than the second 
wall we see no coherent evolution of the domain wall over 
a longer period of time. The motion seen when there is 
no spin density correction present, the LHS of figures [5] 
and 02 is because we have not started the walls in their 
equilibrium position with respect to the anisotropy in the 
system. In this case the domain wall attempts to find its 
lowest energy configuration. What is most clear from 
these figures is that the spin density is able to completely 
distort the profile of the domain wall. Rather than being 
able to shift to its equilibrium configuration the domain 
wall is severely distorted by the spin density torque. Its 
low energy position is now a compromise between not just 
the anisotropy and exchange energies, but also the spin 
density torque exerted on the wall. The low energy profile 
with respect to the spin density contributions would be 



4 




t l f o */ 

FIG. 2. (Color online) A contour plot of the ^-component of the magnetization dynamics M x (t, z) around the second wall 
without an applied current, centered at zo. The left hand figure is without the spin density corrections, the right hand figure 
includes the spin density corrections. The left hand figure shows motion caused by the anisotropy in the wire, on the right 
hand side the spin density waves cause additional distortions, tf = 2.69 x 10 _ s. 
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FIG. 3. (Color online) A contour plot of the z-component of the magnetization dynamics M z (t,z) around the second wall 
without an applied current, centered at Zq. The left hand figure is without the spin density corrections, the right hand figure 
includes the spin density corrections. The left hand figure shows motion caused by the anisotropy in the wire, on the right 
hand side the spin density waves cause additional distortions. One can see that the presence of the spin density readily distorts 
the DW profile. t f = 2.69 x f(r 13 s. 



an oscillating shape, but this is opposed by the exchange 
energy cost it brings. 

Now let us consider what happens when we include the 
current terms in equation The case without any 

correction to the spin density from scattering from a first 
domain wall is shown in the LHS of figure^! here we focus 
on the x-component of the magnetization. One simply 
sees the usual current induced motion of the domain wall. 
The RHS of figure [4] includes the spin density corrections. 
After some time one also sees a distortion of the domain 
wall profile caused by the spin density corrections. The 
new domain wall profile then undergoes coherent motion 
with this new shape. The presence of the spin density 
wave also has a small but noticeable transient effect at 



the beginning of the domain wall's motion, the wall is 
set into motion more quickly than in the case without 
the spin density wave. This effect is short lived and the 
long term motion is, in this case, similar to without the 
spin density oscillations. After a period of time some 
small distortion of the DW profile does appear. However, 
distortions similar to those found when the current terms 
do not contribute are not seen. 

If we lower the applied current density we find the ef- 
fects more pronounced. In fact the presence of the spin 
density appreciably speeds up the motion of the DW 
when a small current is applied, see figure [5J This opens 
up the possibility that in such systems the presence of 
closely packed multiple domain walls assists in the cur- 
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FIG. 4. (Color online) A contour plot of the a;-component of the magnetization dynamics M x (t,z) around the second wall, 
centered at zo- This includes the current terms for the magnetization dynamics, with j e = —2.33 x 10 8 Am ~ 2 . The left hand 
figure is without the spin density corrections, the right hand figure is the full result, showing the accelerated start up to the 
DW motion due to the spin density terms, tf = 6.73 x 10~ 14 s. 



rent induced motion of subsequent DWs, allowing them 
to be shifted with a smaller applied current pulse. The 
spin density wave acts as a local magnetic field, and in 
principle gives one more parameters with which to con- 
trol the current induced domain wall motion. Namely 
the overall and relative amplitudes of the spin density 
correction channels, equation (pTjj) . 

We can also consider what happens for DWs of larger 
widths. For a DW of width L = lOnm the main re- 
sult is still valid. The spin density wave aids the current 
induced DW motion. However if we increase the DW 
width to L = 40nm these effects are already negligible 
for longer timescales. The reason for this is clear, for a 
longer DW the spin density wave is oscillating too much 
over its extension to have a coherent effect. Nonetheless 
for short timescales one still sees quicker initial motion of 
the DW. Furthermore one can see a gradual squeezing of 
the DW caused by the spin density wave, see Fig. [5] The 
edges of the DW are gradually squeezed as they attempt 
to precess in the effective magnetic field of the rapidly 
oscillating spin density wave. This process is limited by 
the increase in exchange energy for a narrower DW. The 
DW ends with a width of approximately 10 — 20nm un- 
dergoing current induced motion. This effect is of course 
absent without the spin density wave. 



IV. DISCUSSION AND CONCLUSIONS 

We have investigated how spin density waves caused by 
spin polarized currents scattering from domain walls af- 
fects the domain wall's magnetization dynamics. In the 
one dimensional limit of atomic chains the spin-torque 
exerted by the spin density waves can be sizable, com- 
parable even with the exchange energy of the bulk fer- 



romagnetic system. For pinned domain walls the inho- 
mogeneous spin distorts the domain wall profile into a 
new shape which minimizes the relevant energies of the 
domain wall. However, when the domain wall is free to 
move through the system the walls are distorted only 
slightly. In this case the main effect of the spin waves 
is to propel the domain wall into motion quicker than in 
their absence. 

The relatively strong effect of the current-induced spin 
density wave on the motion of the domain wall is similar 
to the effect of an external magnetic field. This is be- 
cause the magnetization related to the spin wave appears 
in the equation of motion as an external field, and, on the 
other hand, it is known that by using magnetic field one 
can easily put the domain wall into motion. However, in 
the case of a spin density wave, the effective field acting 
on the domain wall is not homogeneous. On the con- 
trary, it is strongly oscillating. Correspondingly, it can 
affect the short-range interaction of magnetic moments 
making unstable the usual shape of the domain wall. Be- 
ing strongly inhomogeneous, the spin- wave induces some 
displacements and certain disorder of the wall, which is 
important for the DW motion because the first displace- 
ment allows depinning of the wall to put it into motion. 

In this connection, we propose to combine the spin- 
wave-induced dynamics with the short pulses of current. 
It can be realized with a combination of dc current, gen- 
crating the spin density wave, and strong current pulses. 
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FIG. 5. (Color online) A contour plot of the x-component of the magnetization dynamics M x (t,z) around the second wall, 
centered at zo- The left hand figure is without the spin density corrections, the right hand figure is the full result including 
spin density corrections. In this case we apply a smaller current: j e — —5.82 x 10 7 Am ~ 2 . The quicker motion of the domain 
wall due to the spin density correction is now clearly visible in the right hand figure as compared to the left hand figure. 
t f = 6.73 x KT 14 s. 
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FIG. 6. (Color online) A contour plot of the x-component 
of the magnetization dynamics M x (t,z) around the second 
wall, centered at zq. Here we consider a longer DW of width 
L — 40nm. The initial current induced motion of the DW is 
still enhanced by the spin density, for longer times thechange 
in DW velocity is negligible. What is seen is a narrowing of the 
DW width as it moves caused by the spin density wave. The 
applied current is j e = -2.33 x 10 8 Am ~ 2 . t f = 6.73xl(T 14 s. 
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